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Field of the Invention 



The present invention relates generally to receivers of control signals in a communications system, 
such as dual-tone multi-frequency receivers in a telephone network, and, more particularly, to the use of 
5 spectrum estimation techniques in the detection of such signals. 

Background of the Invention 

As a result of tests demonstrating that push-button telephones satisfied consumer desire for faster, 

w easier dialing, the American Telephone and Telegraph Company (AT&T®) began development of its 
TOUCH-TONE signalling system in the early 1950s. Push-button TOUCH-TONE signalling was originally 
devised as a replacement for the once-common rotary-dial pulse signalling system Today, TOUCH-TONE 
signalling controls not only the switching associated with setting up an ordinary telephone call, but also 
many advanced features and functions, such as multiple-party teleconferencing (where, e.g., TOUCH- 

75 TONES are used to control the addition of a party to an on-going teleconference). 

The TOUCH-TONE signalling system uses combinations of "tones" to form control signals representing 
"dialed" digits. The system provides eight tones (or signal frequencies) divided into two groups -- one 
containing four lower pitch tones, and the other containing four higher pitch tones. Each digit of a telephone 
number is represented by two simultaneously generated tones - one tone from each group. Since each 

20 digit is represented by a pair of tones selected from multiple tones, the TOUCH-TONE system is referred to 
as the dual- tone mufti- frequency (DTMF) signalling system. 

The group of lower pitch tones of the DTMF system includes tones of frequencies 697,770,852 and 941 
Hz. The group of higher pitch tones includes tones of frequencies 1209, 1336, 1447 and 1633 Hz. A total of 
sixteen DTMF signals can be produced by combining one tone from each group. For example, the DTMF 

25 signal representing the digit "five" on a conventional telephone keypad comprises simultaneously generated 
tones at 770 Hz and 1336 Hz. Conventional telephone keypads are capable of generating 12 of the 16 
possible DTMF signals (representing digits 0-9 and symbols # and *), while specialized keypads can 
generate all 16 DTMF signals (the conventional 12 plus special symbols A, B, C, and D). To identify a digit, 
a DTMF receiver (or detector) must determine the frequencies present in a DTMF signal and associate with 

30 such frequencies the digit they represent. 

Typically, DTMF signals are not communicated over dedicated, noise-free channels. Rather, they are 
often transmitted over ordinary telephone circuits - -the same -circuits which carry other signals, such as 
voice, music, or data. If these other signals have certain characteristics in common with DTMF signals, a 
DTMF receiver may falsely identify them as DTMF signals. False identification may occur when a signal 

35 carried over a telephone circuit -- for example, a person's voice - -includes one tone (or frequency 
component) from each of the upper and lower DTMF tone groups. 

The problem of falsely identifying voice, music and other signals as DTMF signals is known as talk- 
off. The talk-off problem is not limited to systems which employ DTMF signalling. It affects many types of 
signalling systems which must operate in the presence of noise. 

40 Historically, talk-off has not presented a significant obstacle to the initiation of telephone calls by DTMF 

signalling. This is because, for example, people do not usually speak into their telephone handsets (or 
otherwise intentionally subject telephone circuits to significant noise) when dialing. However, if DTMF 
signalling is to be used in the presence of noise (e.g., voice, as when adding a party to an on-going 
. teleconference), conventional DTMF receivers may exhibit significant rates of talk-off. Thus, talk-off may be 

45 especially troublesome for DTMF signalling occurring during a telephone call, when voice, music or other 
signals are likely to be present. Because of talk-off, many advanced DTMF signal-controlled features and 
functions accessible during a call may be improperly initiated or directed, thereby hindering or preventing 
their beneficial use. 

so Summary of the Invention 



A method and apparatus for "reliably detecting control signals and reducing the frequency of false 
control signal identification, such as that causing talk-off, are provided by the present invention. The 
invention operates to apply tapers (or data windows) to received signals and to determine how welt 
55 transforms of the tapered signals match models of expected control signals. The present invention is 
applicable to the detection of communication system control signals comprising one or more frequency 
components such as, e.g., control signals comprising one or more sinusoidal components (such as DTMF 
signals, call progress signals (e.g., ring and busy signals and dial tones), supervisory audio tone signals 



2 



EP 0 550 144 A1 



(used in, e.g., cellular telephony), and the CCITT (International Telephone and Telegraph Consultative 

Committee) No. 5 signalling system, among others). 

In an illustrative embodiment of the present invention, digital signals which may contain DTMF signals 

to be detected, as well as other signals, are received and sliced into time segment s which may be 
5 o yeria P.Qiflg. To each of a plurality of copies of a segment is applied one of a plurality of tapers, such as 

one of the Discrete Prolate Spheroidal Sequences (DPSS). Discrete Fourier Transforms of these tapered 

segments are determined and combined to form one or more estimates of the distribution of signal energy 

as a function of frequency. .A ratio of best-fit DTMF sinusoidal model energy to model misfit energy is 

formed for each of the DTMF sinusoidal frequencies. If the value of the ratio exceeds a threshold, there is 
io likely to be a sinusoid at the DTMF tone frequency associated with the quotient and not elsewhere. If 

sinusoids are found at certain pairs of frequencies (and not elsewhere), DTMF signals may be identified. 

Additional tests concerning, e.g., tone duration and spacing, may be included in the DTMF signal 

identification procedure. 



75 Brief Description of the Drawings 



Figure 1 presents an illustrative telecommunication system comprising a plurality of individual telephone 
networks. 

( Figure 2 presents an illustrative control function for the telephone network of Figure 1. 

^ 20 Figure 3 presents a portion of the control function presented in Figure 2 concerned with receiving and 

detecting control signals. 

v Figure 4 presents the eight lowest-order discrete prolate spheroidal sequences determined for NW = 4. 

Figure 5 presents the eigenvalues of discrete prolate spheroidal sequences determined for NW = 4. 
j Figure 6 presents an illustrative DTMF receiver according to the present invention. 
25 Figure 7 presents a time progression of signal segments provided by an illustrative embodiment of the 

present invention: 

•* Figure 8 presents the embodiment presented in Figure 6 illustratively augmented to provide a maximum 
frequency variation check, a maximum amplitude variation check, a maximum phase variation check, and a 
twist check. 

30 Figure 9 presents an illustrative embodiment of the present invention which performs a simultaneous 

estimation of signal amplitude and frequency. 

Figure 10 presents an illustrative embodiment of the present invention which provides for the simulta- 
neous estimation of the amplitudes of two component tones of a DTMF signal. 



35 Detailed Description 

A. Introduction to the Illustrative Embodiments 

A communication system comprises a multiplicity of devices or systems which need to communicate 

40 with each other, and a set of one or more channels, coupled to the devices or systems, which facilitate 
communication. Perhaps the most well-known and widespread of all communication systems are the 
telecommunication systems, such as telephone networks. These include, for example, local- and long- 
distance telephone networks (e.g., the AT&T long-distance telephone network), and private branch ex- 
change networks .(PBXs) (e.g., the AT&T DEFINITY® PBXs) used to provide special services to many 

45 businesses and institutions. A communication system device which transmits or receives information over 
the system is referred to as a system station. 

Figure 1 presents an illustrative communication system comprising a plurality of individual telephone 
networks and a collection of stations. The individual telephone networks include PBXs 3, providing features 
such as a central answering position, station-to-station intercom calling, and connectivity to local and long 

so distance telephone networks; a local-telephone network 5, providing not only services associated with 
completing local telephone calls and connecting local stations to a long-distance network 6, but also integral 
services such as, voice-messaging, three-way calling, call-waiting, etc.; and a long-distance network 6, 
linking individual local networks 5 together and also providing integral services. Stations may include 
telephones 1, facsimile machines 2, and personal computers 4 (PCs), among other things. The illustrative 

55 telecommunication system of Figure 1 may further comprise, e.g., network switching systems (e.g., the 
AT&T 4ESS™ and 5ESS® switches) network-based teleconferencing (e.g., AT&T Alliance® Teleconferen- 
cing) and messaging (e.g., AT&T Message Service), PBX-based messaging and switching systems, voice 
transaction systems, etc. Components (or parts thereof) of a communication system may be referred to 
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generically as system elements. 

In order for one system element to communicate with another, a system element referred to as a 
channel must be established. A channel may take the form of a dedicated coupling between elements or a 
shared coupling used by several pairs or groups of elements simultaneously. A channel may comprise an 

5 electrical coupling, a radio link, a fiber-optic element, an acoustic coupling etc 

To establish and maintain communication channels (among other things), communication systems (such 
as telephone networks) may comprise one or more elements referred to as system control functions A 
system control function provides control of a communication system element in response to control sionals 
communicated over a channel by another system element. For example, a system control function may 

w control system sw.tch.ng capability directing information flow through a telephone network That is a 
communication channel may be established when a station desiring to initiate communication (a source 
station) with another station (a target station) transmits control signals over the network identifying the 
target to one or more system control functions. A system control function, in response, directs the coupling 
(or sw.tchmg) of channels to facilitate the desired communication between source and target stations 

15 <™ ln i a f dd,, '° n t0 ' 0r instead of - controlling a switching capability, a system control function may control 
special functions, services, devices, and systems ("systems") which augment or enhance a communication 
system s basic switching and communication functions. Like a switching capability, these special systems 
are controlled by a control function in response to control signals. In telephone networks, these special 

20 ZsZt m T,' n « t Ca "- f0rwardin9 and - trans,er s ^™- voice-messaging, teleconferencing, transaction 
20 systems, etc. (,n effect, any function, service, device or system which may be initiated or directed by 
control signals). ' 

Figure 2 presents an illustrative system control function for the telecommunication system of Figure 1 
This control function" could be resident anywhere in the telecommunication system where control of 
switching or special systems is needed (e.g.,: at a station, at a PBX, in a channel, in a local telephone 
25 network, or ,n a long-distance network). Input to the system control function (CF) are a plurality of network 
communication channels 7, from a plurality of system elements (not shown). These input channels are also 
coupled to, e.g., a network switch 8, Switch 8, in turn, is coupled to other (output) channels 9 leading to 
other system elements (also not shown). Switch 8 may perform the coupling of any. of the input channels 7 

,» n nMh°J ° T T S ( • maV Pr ° Vide Sp6Cial funCti0ns ' a " under the direction of control Processor 
30 10 of the system control function. 

Control processor 10 of the system control function receives control signals from a system element via 
input channel 7 and mput control signal interface 13. Control signals may comprise, e.g., electrical or optical 
signals representative of a telephone number of a target station or a code for a special function desired by 

,< I T I ? S ' 9na • WhGn detGCted bV COntr °' Si9nal interface 13 ' are P assed to control processor 10 
35 to be used as input mformation for software executed by the processor 10. This software, stored in program 
memory 11, provides the functionality associated with the system control function. As a result of executing 
the software,, control processor 10 determines other signals necessary to, e.g., facilitate connection of 
source and target stations, or to provide desired special functions. (To perform its tasks, control processor 
10 may emp oy a scratch pad memory 12 to record and accumulate information during processing It may 
40 a so prov,de the network with control signals, via output signal interface 14, to further facilitate establishment 
or a communication channel or to provide special functions.) 

rn JrniV° ? un f c rst00d that tne number - and configuration of hardware components for the system 

control function of Figure 2 is merely illustrative. Other component configurations are possible. For example 
a control func ,on not comprising a central processor unit may be provided. Naturally, the present invention 
is applicable to any control function configuration which requires the detection of communication system 

control sicjnsis. 

Control signals received by a telecommunication system control function may take the form of DTMF 
Signals. As discussed above, DTMF signals may be communicated over telephone circuits and therefore 
may be ,nc uded with other "noise" signa.s, such as voice, music, computer data signals, physical noise 
and crosstalk from adjacent circuits. The input control signal interface (l/F) 13 is charged with the task of 

tZT?T, T T tf01 Si9nalS Pr6SenCe ° f thiS n0iSe ' ln P erf °™ng this task, the control signal 

.nterface 13 should not mistake such noise for valid DTMF control signals (talk-off). To detect true DTMF 
control s.gnals and avoid talk-off, an input control signal interface may include an illustrative DTMF signal 
detector according to the present invention, as shown schematically in Figure 3 

55 Figure 3 presents a portion of the control function presented in Figure' 2. Signals, x(t), from input 

channels 7 are provided to input control signal interface 13 and to an illustrative DTMF signal receiver 13' 

• ^cording to the present invention (other circuitry and components of interface 13 are excluded from Figure 
3 for clarity). As discussed above, signals. x(t). may be. e.g., physical electrical or optical communication 
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signals, typically in digital form, which may include DTMF signal information intended to inspire actual 
system control, as well as noise (if x(t) is in analog form, it is common that the signals will be converted to 
digital form prior to DTMF detection). It is the task of illustrative DTMF signal receiver 13* to receive signals, 
x(t), and to correctly isolate or identify any actual DTMF control signals which may be present in x(t). As a 
• 5 result of identifying a DTMF signal in x(t), illustrative DTMF signal receiver 13' produces other electrical 
signals, g(t),. which represent identified DTMF signals present in x(t). Signals g(t) are then provided to 
control processor 10 to inspire desired network control. (The switch and special functions described above 
for telephone, networks are illustrative of systems or devices controlled by application of physical control 
signals. In light of this, the ordinary artisan will understand that the present invention is applicable to many 
w systems and. devices which require the identification of physical control signals in the presence of noise.) 

For clarity of explanation, the illustrative embodiments of the present invention are presented as 
comprising individual functional blocks (including functional blocks labeled as "processors'). The functions 
these blocks represent may be provided through the use of either shared or dedicated hardware, including, 
but not limited to, hardware capable of executing software. (Use of the term "processor" should not be 
75 construed to refer exclusively to hardware capable of executing software.) Illustrative embodiments may 
comprise digital signal processor (DSP) hardware, such as the AT&T DSP16 or DSP32C, and software 
performing the operations discussed below. Very large scale integration (VLSI) hardware embodiments of 
the present invention, as well as hybrid DSPA/LSI embodiments, may also be provided. 
) As part of the process of identifying DTMF signals, the illustrative embodiments of the present invention 

20 estimate energy distributions of received digitized control signals. A signal's energy spectrum represents 
the distribution of signal energy as a function of frequency. In general, when determining energy spectra of 
received control signals, it is desirable to do so in a way which minimizes any mismapping or spurious 
leakage of . energy from one frequency in the spectrum to another. If such leakage is unchecked, a 
spectrum estimate may misrepresent the energy distribution of a signal and, consequently, obscure true 
. 25 spectral features. To limit leakage in spectrum estimates, the illustrative embodiments of the present 
j invention employ data windows or tapers. Tapers have the effect of concentrating energy around 

t appropriate frequencies in a spectrum (i.e., where the energy actually is), thereby reducing the distortion of 

true spectral features. 

30 1. Discrete Prolate Spheroidal Sequences as Tapers 

Several of the illustrative embodiments of the present invention employ special tapers referred to as 
Discrete Prolate Spheroidal Sequences (DPSS). These DPSS are especially good at reducing energy 
leakage in spectrum estimation. While other types of tapers may he employed in embodiments of the 

35 present invention, because of the energy-concentrating properties of DPSS tapers, it may be instructive to 
consider them in some detail. 

To simplify notation in this and the following sections, time t has been designated as an index taking on 
integer values and a unit sampling interval has been assumed, so that t = 0, 1, 2,..., N-1. The frequencies in 
the formulae which follow, such as fo, are therefore numbers in the Nyquist band (-1/2, 1/2]. The bandwidth 

40 W of a frequency band is therefore a fraction of the width of the Nyquist band, However, the 

physical frequencies and bandwidths of the embodiments discussed below are given in Hz, as is 
conventionally done. To convert a physical frequency of 200 Hz to a frequency fo as appearing in the 
formulae below (where data is sampled at a rate of 8000 Hz), f 0 is calculated by dividing the physical 
frequency 200 Hz by the sampling rate, 8000 Hz: 

45 

200 Hz 

fo ~ 



8000 Hz 



50 



Similarly, a bandwidth of 1000 Hz corresponding to a bandwidth W appearing in the formulae below is 
calculated as: 



55 = 1000 Hz 



8000 Hz 
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Consider a finite duration sinusoidal signal 



A 0 e 



having a complex amplitude, A 0l and a frequency, f 0t sampled at discrete times t = 0,1 ,2,...,N-1 . A Discrete 
Fourier Transform (DFT) of the signal multiplied by a taper, w(t), is: 



10 



v(/) = N ± X A 0 e- t2nfot w(t) , 

: = 0 



(.4.1.1) 



and the energy of the tapered sinusoid is \y(f)f. The taper, w(t), may be chosen to concentrate most of the 
75 energy of the sinusoid within a restricted frequency band, (f Q -W t f 0 + IV], of a spectral estimate. However, 
some of the energy of the sinusoid may leak to other portions of the Nyquist band. The total energy of the 
sinusoid within the Nyquist band, E r , is: 



20 



E r = J \y(f) \ 2 df , 

-1/2 1 1 



(A. 1.2) 



25 and is the same as the total energy of the tapered signal, 



30 



35 



V Ijc(0w(/)I 2 . 

1 = 0 

Similarly, the energy in the restricted band (f 0 -W, f 0 + W], E R , is: 

; 



E * = J \y(f)\ 2 df. 



M:L3) 



The ratio of energy within the restricted band to the total energy in the Nyquist band is referred to as 
40 the concentration ratio, X: 



45 



fo-W 



df 



1/2 

/ 

-1/2 



y(f) 



(AAA) 



df 



50 



55 



where X, quite naturally, has a maximum theoretical value of one. The better a taper is at concentrating 
energy in the restricted band, (i.e., reducing leakage), the closer x is to one. Because reduction of energy 
leakage is desirable for accurate spectrum estimates, "good" tapers are those which maximize X. The 
tapers maximizing equation (A.t.4) are several which satisfy the symmetric, Toeptitz eigenvalue problem: 
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N ' 1 sin iKW(t-t') 



z 



w(r') = Xw(t) , 



(4.1.5) 



where t = 0,1, 2.....A/ - 1. The solutions to this problem are the DPSS, v k (t), /c = 0,1,2 A/ - 1. 

A concentration ratio \ k is the eigenvalue corresponding to the eigenvector of order k, v k {t). The order of 
the DPSS v k (t) is defined so that 1>\ 0 >\i>\ 2 ... X w _ ^0. Moreover, the [2A/W J lowest-order DPSS have 
eigenvalues close to one, and therefore provide good energy concentration within the restricted band 

70 (where |_*J indicates the integer portion of *). A plurality of the DPSS, » k (t), typically O^k^K - 1 where 
Q<K<>\_2NW\ t may be used in combination to provide an accurate estimate of the spectrum of a signal. 

Figure 4 presents the eight lowest-order DPSS determined for NW = 4. Figure 5 presents eigenvalues of 
these DPSS as a function of k. Figure 5 shows that for k£7, the DPSS have eigenvalues, X*, which are 
greater -than or approximately equal to 0.7. Moreover, for k£5, \ k ^^. As such, these first eight DPSS 

is generally, and the six lowest-order DPSS specifically, are very good at concentrating energy and preventing 
leakage in spectrum estimates. Not all of the DPSS are this good, as Figure 5 shows for K>Q (and 
especially for k^10). However, for ks\_2NW\ -1. values of \ k are generally acceptable for purposes of 
energy concentration in spectrum estimation. 

20 2. Numerical Evaluation of the DPSS 

The DPSS obey the following symmetric tridiagonal eigenvalue equation: 



25 



£ S tu v k (u) = e*v*(0 , 

u=0 



(A. 2.1) 



30 



where: 



35 



O 



' tu 



40 



(t 0 -t) 2 co$2nW 
y(r + l)(N-l-0 



w=f-l 

u-t . 
u = t + 1 
U-rl>l 



(A. 2. 2) 



for t,u = 0,1,2,..., N-1, where 



45 



jV-1 



and d k is another eigenvalue associated with v k (t) t as discussed by D. Slepian, Prolate Spheroidal wave 
functions, Fourier analysis and uncertainty - V: The discrete case, 57 Bell System Tech. J., 1371-1429 
so (1978). The DPSS can be found for a wide range of values of W and N by solving the pair of 



55 



N+l 



tridiagonal matrix eigenvalue equations of the form (A. 2.1) using conventional procedures such as that 
provided by Smith et al., Matrix eigensystem routines - EISPACK guide, Lecture Notes in Computer 
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10 



15 



Science No. 6 (the complete NxN system need not be solved because the DPSSs are either odd or even 
about their midpoints at t= t 0 ). Because the DPSS are scaled to have unit energy, 

Z (v*(f)) 2 = l, 

the concentration ratios \ k , k = 0,1,2 N-1, can be determined from the sum: 



= L v k (t) \ ' v k (t ) , (A. 2.3) 

M '=0 K(t —t ) 



3. Tapers Other Than the DPSS 

In some applications it may be desirable to use sets of tapers other than the DPSS. For example, there 
20 are alternative sets of tapers with energy concentration properties which are almost as good as those of the 
DPSS, such as truncated products of Hermite polynomials with a Gaussian density. 

A second class of alternative tapers may be obtained by replacing (A.1.3) in the discussion above with: 

E R = J |y(/)| 2 5(/) df, 

-Yi 



where the weighting spectrum S(f) might, for example, contain the average out-of-band energy distribution. 
30 This may provide enhanced protection against sources of talk-off having this energy distribution, and may 

be preferred for environments with identifiable noise characteristics. 

A third class of alternative tapers is given by the discrete Karhunen-Ldeve expansion. If process 

statistics are known accurately, tapers given by the Karhunen-Ldeve (K-L) expansion have the fastest 

convergence, in a statistical sense, of any possible taper set and give maximum-likelihood estimates. 
35 If the Karhunen-Loeve expansion is converted to the frequency domain and the frequency domain is 

restricted as in (A. 1.3), the tapers are solutions of the integral equation: 



JO 

where S(£) is the power spectrum of the signal causing talk-off. (If S(£) is constant in (/ 0 - W t f 0 + W], the 
K-L tapers are frequency shifted DPSS.) Linear combinations of the discrete prolate spheroidal sequences 
45 efficiently approximate the restricted K-L tapers. The hybrid hardware described earlier may be adapted to 
use this statistical information, wherein the eigen-coefficients are computed by VLSI circuitry and the 
subsequent processing performed by a DSP. This approach has the advantage that it may be made 
adaptive: using either the DPSS or the band-limited Karhunen-Loeve tapers for 'average' speech initially, 
the eigencoefficients could be transformed to reflect an individual's speech statistics. 

50 

B. An Illustrative Multiple Taper DTMF Receiver 

An illustrative DTMF receiver 13' according to the present invention for use in, e.g., an input control 
signal interface 13, is presented in Figure 6. The receiver comprises a signal slicer 15 (which receives input 
55 signals which may contain DTMF signals and provides usable segments thereof), a DFT processor 20 
(which applies tapers to the signal segments and computes the DFT of the tapered segments), a memory 
25, coupled to DFT processor 20 (for storing taper data used by DFT processor 20), a spectrum amplitude 
estimator 30 (for estimating the amplitude of a spectrum at specific frequencies), a statistical F-Test 

8 
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processor 40 (for determining the how well the energy of the received signal is explained by a sinusoidal 
model at a given frequency), and a DTMF signal identifier 50 (which, depending on information received 
from F-Test processor 40, may provide an output indicative of a particular DTMF signal to a system or 
device 60 which operates on received DTMF signals). Illustrative software related to this embodiment 

5 (written in the "C" programming language) is presented in an Appendix hereto. This software may be 
executed on conventional hardware platforms, such as, e.g., a SUN® SPARCstation 1.® 

Signal x(t) is a discrete series of values which may represent, among other things, a DTMF signal. The 
signals is input to signal slicer 15 which divides this discrete series into sets of contiguous or^verjappjng 
signal segments, each comprising N elements. That is, x* s) (t), $ = 0,1,2,.... Figure 7 presents a time 

70 progression of the segments, x^ s) (t) f for signal x(t). A signal segment may be contaminated with noise, and a 
DTMF signal may begin or end any time within a segment. To detect short duration DTMF signals and 
short -inter-digit intervals, a signal segment may be chosen to be less than 30msec in duration; an 
illustrative value is 10msec, corresponding to a segment of about N = 80 samples at the standard 8000Hz 
telephone system sampling rate. 

/5 DFT processor 20 receives signal segments from signal slicer 15 and computes 24 DFTs of length N 

from each segment, x< s) (f)- Each DFT is of the form: 

A s) (fj) = N £x M (i)VkU)e- i2nf >\ (5.0.1) 

v nr*A-X *A "TIT /-<k->v ■ bu*-*-*^ ^ - 
where f h 1 ^ j £ 8, represents the eight DTMF tone frequencies, and v k (f), 0 ^ k £ 2, represents the three 
lowest-order DPSS specified by values for N and W. The resulting values 



20 



25 



yi'Hfj) 



are referred to as eigen- coefficients and are provided as output from DFT processor 20 to spectrum 

30 amplitude estimator 30. (DFT processor 20 may be replaced by a processor performing other transforms 
such as Laplace, Walsh, etc.) The values of the tapers y k (t) may be stored in memory 25, which may be, 
e.g., a semiconductor programmable read-only memory (PROM), and recalled as needed. There is thus a 
DFT for each of the eight frequencies at which DTMF tones are found and each of the three lowest-order 
DPSS, yielding 24 DFTs for each segment x< s, (0- 

35 The number of tapers used by the DFT processor 20 is chosen to be the integer portion of 2NW - 1. 

The tapers are determined based on a set of frequencies that include a spectral feature of noise ("noise" 
being non-DTMF signals) in addition to that which contributes to talk-off. Thus, the bandwidth 2W is chosen 
wide enough to discriminate between DTMF signals and other signals which may cause talk-off. In this; 
embodiment, a value for W should be chosen such that there can be at most one frequency component 

40 from a valid' DTMF signal in the range of frequencies (f, + W, f } - W] for 1 S j £ 8. A different value for W 
may be chosen for each f h An illustrative value for W for all fj is 230 Hz/8000 Hz. This allows use of the 
three lowest-order DPSS (the integer portion of 2NW - 1 equals 3, where N = 80) and defines bandwidths 
which avoid dial tone frequencies used in most of North America (for other values of W, it may be useful to 
eliminate dial tone frequencies with a high-pass filter. 

45 To reduce the number of computations required to determine 



so equation (B.0.1) can be implemented by successively summing precomputed values retrieved from memory 
25. These precomputed values comprise the application of tapers to all possible values of data. For 
example, assuming x(t) is an 8-bit u-law signal, x(t) may take on any of 2 8 (or 256) possible values. 
Memory tables for the expression 



55 
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may be precompiled and stored in memory 25 for each k, t, and frequency f h Values for j, k, t and x(t) may 
be used to index the tables and retrieve the precomputed values. Since there are only eight DTMF tone 
frequencies and, for example, three values fork required, only 8x3x2xA/ = 4SN tables of 256 entries need 
to be stored (where N is the length of the segment and the factor of two arises due to real and imaginary 
5 components). Summing precomputed values to form 

yi'Hfj) 

10 will alleviate much of the computational burden of performing real-time DFTs. . 
Amplitude estimator 30 receives the eigen-coefficients, 



75 



20 



and determines an estimate of the complex amplitude of a transformed signal segment at the eight defined 
frequencies at which DTMF tones may be found. At each such frequency, f h an amplitude estimate in the 
s ih segment, A {s) , is formed as a weighted sum of the eigen-coefficients: 



/ l'^(0)yit J) (//) 



25 



n ■ 



k = 0 



(B.0.2) 



where 



30 



Vkif) = N £v k (t)e- i2n f'\V k (0) 
r =0 



35 is the DFT of » k (t} evaluated at f = 0) and K = 3. 

Expression (B.0.2) for estimated amplitude, A {s) (fj), is derived by minimizing an expression for model 
misfit energy: 



40 



0 K ' 1 



(5.0.3) 



45 



with respect to A at each frequency / = f { . Model misfit energy in the frequency domain represents how a 
transform of an actual signal, y*(//), differs from a transform of a sinusoidal model of the signal, A V k {0), 
where 



50 



V k (0) = V v*<0. 

t =0 



55 



If the signal is accurately represented by a sinusoid at f h e 2 (fj) will be small. If not, e 2 (f,) will be large. The 
sinusoidal model is useful since DTMF signals to be detected comprise sinusoids. A(f) is determined by 
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de 2 (fj) 



0, 



10 



and solving for A. 

The values of A (s) (/y), 1 ^ j < 8, are provided as output from the amplitude estimator 30 to the statistical 
F-Test processor 40. 

The statistical F-Test processor 40 determines which of the eight DTMF tone frequencies, if any, are 
presentrin a group of R successive signal segments in question. It does this by determining values for a 
test variable, F s (f,-): 



15 



20 



25 



30 



35 



o 



40 



45 



50 



* (s) 

A (fj) 



2 K- \ | 

Z Kco) 

k=0 ' 



F s(fj) = 5~</?-l) K-l ~ TT 

X Z yi s) (fj) - a { \fj)v k {0)\ 

s=S k=0 1 



(5.0.4) 



where R = 3, K = 3 and S is the index of the latest segment. The numerator of the expression for F s (/,) 
represents the energy of a best-fitting sinusoidal model at frequency f h while the denominator represents 
the energy associated with the difference between the received signal and the sinusoidal model, sometimes 
referred to as the misfit energy. If the group of signal segments is well-approximated by the model (i.e., if 
the received signal actually contains a signal sinusoid at frequency fj), the numerator of F s (fj) will be nearly 
the signal energy, the misfit energy will be small, and the ratio, F s (/,), will be large. Conversely, if the model 
does not match the signal, the modeled (or explained) energy will be small, the misfit energy will be large, 
and the ratio F s (fj) will be small. 

The F-Test discussed above (B.0.4) is highly frequency sensitive. That is, F s (ff) will tend to be large 
only when a group of R segments includes a sinusoid at a frequency very close to an expected frequency, 
Because of DTMF tone generator component tolerances, aging, etc., and. telecommunication system- 
based frequency shifts, DTMF signals may not be communicated exactly at the expected frequencies, f h 
Consequently, reduced F-Test frequency sensitivity may be desirable. Moreover, it may be needed to 
satisfy DTMF system specifications which may require a DTMF receiver to identify tones which are within a 
certain frequency range about expected frequencies. A segment length of 10 msec provides a relaxation in 
the F-Test frequency sensitivity, when compared with longer length segments, to allow for variation of ±- 
(1.5% +2Hz) of expected DTMF frequencies, f h 

F-Test processor 40 determines values for F s (f,) at each of the four low DTMF tone frequencies (697 
Hz, 770 Hz, 852 Hz, and 941 Hz) and each of the four high DTMF tone frequencies (1209 Hz, 1336 Hz, 
1447 Hz, and 1633 Hz). Values of F s (fj) are output from F-Test processor 40 to a DTMF tone identifier 50. 

The DTMF tone identifier 50 receives the values F S (A / ) and determines whether such values indicate the 
presence of a DTMF signal. A DTMF signal may be identified when one and only one determined value of 
F s (fj), 1^y^ 4, and one and only one determined value of F s (/,*), 5 £ / £ 8, exceed a selected confidence 
threshold. 

In the absence of a valid DTMF signal, the random variable F s (fj) is approximately distributed as a 
central "F" distribution with 2R and 2R(K - 1) degrees of freedom. The confidence thresholds for values of 
F corresponding to a given number of degrees of freedom can either be calculated or found in a set of 
statistical tables (e.g., Abramowitz and Stegun, Handbook of Mathematical Function With Formulas, 
Graphs and Mathematical Tables, Ch. 26(1964)). The DTMF tone identifier 50 can therefore determine, 
with a desired level of confidence, whether the energy present in a spectrum can be attributed to a sinusoid 
at a given frequency. For example, if K = 3, R = 3, and F s (//)>4.82, a confidence level of 99% is indicated. 

If one and only one value of F s (fj) associated with each tone group exceeds a confidence threshold, a 
DTMF signal is identified by DTMF signal identifier 50 as the ordered pair of frequencies, 



55 



(fj L Jj H \ 



at which F s (f,) exceeds the threshold (where 
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is the one frequency at which F s (/,), exceeds the confidence threshold, and 




is the one frequency at which F s (f,), 5^8, exceeds the same threshold). Responsive to identifying a DTMF 
signal, an indicator signal specifying the DTMF signal, g m , 0^/77^15, is passed to system element 60 which 
operates in response to the detected DTMF signal. In the illustrative control function presented in Figure 2, 
system element 60 comprises control processor 10. 

In generating signal g m , DTMF tone identifier 50 may enforce certain requirements regarding DTMF 
signal length, DTMF signals in consecutive segments, and time gaps between DTMF signals defined by a 
DTMF system specification. For example, a DTMF system specification may require the rejection of any 
DTMF signal which is not at least 40 msec in duration. A DTMF system specification may require that when 
a given DTMF signal is present in consecutive, segments, g m should reflect a single DTMF signal and not a 
DTMF signal Per segment. Also, a system specification may require that when time between consecutive 
like-DTMF signals is less than a threshold (e.g., 20 msec), the consecutive like-DTMF signals should be 
treated as a single DTMF signal. 

The specific characteristics and form of g m may be determined by the logical and electrical interface 
requirements of system element 60. In the context of an illustrative network control function discussed 
above, these interface requirements may be dictated by control processor 10. 

If other than one and only one value of F s (f,-) for each tone group exceeds the threshold, then DTMF 
tone identifier 50 will not pass any indicator signal to system element 60. 

Regardless of whether a DTMF signal is detected in one or more signal segments, the entire detection 
process continues in real time. That is, the operation of signal slicer 15, DFT processor 20/memory 25, 
amplitude estimator 30, F-Test processor 40, and DTMF signal identifier 50 is continued to determine 
whether subsequent signal segments, x< s) (f). contain DTMF signals. To aid multiplexed applications, DTMF 
tone identifier 50 may provide a signal, NS, to signal slicer 15 to provide control of the application of signal 
segments, x< s) (0, to DFT processor 20 for processing. 

Although each calculation of F s (fj) by F-Test processor 40 requires DFT and amplitude information for 
each of R consecutive segments, values for 

yP 

and A (s) , once determined by DFT processor 20 and amplitude estimator 30, respectively, may be stored by 
F-Test processor 40 for use in subsequent calculations of F s (fj)~. 

The illustrative embodiment discussed above and presented in Figure 6 may be augmented to provide 
additional testing of the received signal segments. Such additional testing may be desirable to determine 
the presence of DTMF signals with even greater certainty, or to satisfy criteria for "good" DTMF signals 
which may be defined by a DTMF system specification. Figure 8 presents the embodiment of Figure 6 
illustratively augmented to provide a maximum frequency variation check 51, a maximum amplitude 
variation check 52, a maximum phase variation check 54, and a twist check 56. This embodiment also 
includes an interface 58 to provide the requisite logical and electrical compatibility with system element 60, 
as well as the processing associated with DTMF signal duration, consecutive segments with like-DTMF 
signals, and inter-signal gaps. ' 

1. Maximum Frequency Variation of DTMF Tones 

A DTMF system specification may place requirements on the degree to which DTMF tones can vary 
from agreed-upon standards and still be accepted as valid. For example, a system specification may 
provide that a received sinusoid with frequency f shall not be accepted as a valid DTMF tone (with 
frequency if \fj - ^0.035 f h This requirement can be accomplished by making amplitude estimates not 
only at frequencies f h but also at frequencies 1.035f, and 0.965//. 

As before, a check is made by DTMF tone identifier 50 to see if there is one and only one tone in each 
of the low and high tone groups. If not, DTMF tone identifier 50 may reject the signal as a non-DTMF signal 
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and notify signal slicer 15 as described above. If there is only one tone in each of these groups as 
determined by DTMF tone identifier 50 (designated as 



ft a™*/" 



for frequencies in the low and high tone groups, respectively), a check is made by processor 51 to 
determine whether the estimated amplitudes at both 

• //and// 

are greater than the estimated amplitudes evaluated at frequencies offset by i 3.5% of 

//and//, 

respectively. If so, actual received tones at f L and /h are within the ± 3.5% requirement and, as such, may 
be accepted as constituting a valid DTMF signal at 



75 



20 



fL f H 



25 



If not, then the requirement is unsatisfied and no DTMF signal at 

tL rH 



J J *Jj 



may be recognized. Signal slicer 15 may then be notified by process 51 (via the NS signal) to provide the 
30 next signal segment This test may be conducted for each of the R segments considered (see Appendix). 

2. Maximum Amplitude Variation of DTMF Tones 

DTMF signal tones exhibit relatively constant amplitude. To help insure the detection of actual DTMF 
35 tones, an amplitude consistency check for the R segments considered by F-Test processor 40 may be 
performed illustratively as follows. 

The above-described embodiment of the present invention determines an amplitude estimate for each 
(^) segment at each DTMF tone frequency, ^ <s> (/,*). Assuming that one tone from the low and high tone groups 

with frequencies 

40 



ft and //, 



respectively, have been identified by DTMF tone identifier 50, an average absolute amplitude over R (e.g., 
45 3) segments at each frequency is 

// and// 

so determined as: 



1 S~{R-\) 



55 S = S 



The value 
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w 



A(fj LH ) 



is used by processor 52 to calculate at each frequency a quantity 

FAMPiff""). 

which compares average amplitude to an amplitude dispersion between segments: 

(R-lKA(fj LH )) 2 



FAMPUt'") = 



)5 



S-(R-l) 



T - 



(5.2.2) 



Because 



20 



P AMP(fi LM ) 



is inversely proportional to the dispersion of 



25 



30 



35 



40 



small values for 

FAMPiff") 

are indicative of a tone with significantly varying amplitude over R segments. Tones 

which vary beyond the degree expected for DTMF tones, but which might otherwise satisfy the F-Test, are 
rejected by processor 52. Thus, if either of the tones at 



fj orfj 



fail this check, no DTMF signal at 



45 



L f H 



may be identified. Under such circumstances, signal slicer 15 may be notified by processor 52 to provide 
50 the next signal segment- A threshold value for 

FAMPiff"") 

55 may be chosen which, if not exceeded, will require rejection of the DTMF signal. An illustrative value for the 

FampUj LM ) 
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threshold is 2.5. 



3. Maximum Phase Variation of DTMF Tones 



10 



As with tone amplitude, the phase of a DTMF tone is relatively constant for its duration. To help insure 
the detection of actual DTMF tones, a maximum phase variation check for the R segments may be 
performed by processor 54 illustratively as follows. 

A value for D (where R = 3) is determined by processor 54 for each of the tones identified by F-Test 
processor 40, 



//• and//, 



75 



20 



such that 



£>(//■") = 



) 



A (TT (ff" ti ) 



where values for 



25 



are complex. (For R>3, similar formulae may be used.) The greater the phase shifts which may be present, 
the greater the value of D. Tones exhibiting phase shifts which vary beyond the degree expected for DTMF 
30 tones, but which might otherwise satisfy the F-Test, are rejected by processor 54. If either of the tones at 



ft °rf? 



35 fail this check, no DTMF signal at 



o 



f L f H 



40 may be identified. As described above, signal slicer 15 may be notified. A threshold value for 



45 may be chosen which if exceeded shall require rejection of the DTMF signal. Illustrative values for the 
threshold for D are 0.33 and 1 .0. 

4. Minimum and Maximum Twist of a DTMF Signal 

so The specification of a DTMF signalling system may set fairly wide limits on the permissible relative 

amplitudes of the two tones in a received DTMF signal. These wide limits allow for a phenomenon known 
as twist, caused by imperfect tone generators in telephone sets and frequency-dependent attenuation in a 
communications network over which tones are transmitted. If A L is the amplitude of a. lower frequency tone, 



55 




identified by tone identifier 50 and A H is the amplitude of a higher frequency tone, 
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w 



15 



20 



45 



50 



identified by tone identifier 50, a DTMF signalling specification may require that 



-6dB < -r—^, < SdB, 
\A H \ 



if the pair of tones 



f.L f H 

is to be identified as a DTMF signal. Twist check processor 56 determines whether the value of 

\A L \ 
\Ah\ 



falls outside the specified bounds. If so, then the DTMF signal associated with these tones may be rejected 
25 and signal slicer 15 may be notified 

C. An Illustrative Embodiment Providing Simultaneous Estimation of Amplitude and Frequency 

As discussed above, convention provides that an ideal DTMF signal comprises a pair of sinusoids at 
30 defined frequencies. As a practical matter, however, tones, which form real DTMF signals may deviate 
slightly from the definitions. Such deviations may be caused by, among other things, imperfectly-tuned 
DTMF signal generators. Naturally, specifications for DTMF systems may require. that small deviations from 
defined frequencies be accommodated by DTMF receivers. In a embodiment discussed above, small 
deviations are accommodated through use of short signal segments which have the effect of lessening the 
35 frequency selectivity of the F-Test. Below, a further embodiment for accommodating deviations from 
defined frequencies is described. This embodiment operates to determine an estimate of frequency 
deviation, in addition to tone amplitude, and reject tones which deviate in frequency beyond an acceptable 
amount. ' 

Figure 9 presents an illustrative embodiment for the contemporaneous estimation of amplitude and 
40 frequency. Signal slicer 15, DFT processor 20, and memory 25 are the same as those discussed above. 
Amplitude and frequency estimator 70 receives eigen-coefficients 

. >i s) (fj) 

from DFT processor 20 and determines A (s) (fj) as per equation (B.0.2). Using the estimated A< s) (//). estimator 
70 determines an estimate of frequency deviation about 

as follows: 



55 
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a ; = 



*j:G k yi'\fj) 



Jt =0 



' * = 0 



(C.0.1) 



where 



/o 



/V - 1 



N - 1 



G* = -i2tc X (t-t 0 )v k (t) t H k = -2ti z X (f-f 0 ) v A (0. 



75 and 



20 



25 



N-l 

Alternatively, an estimate of the frequency offset a/ at segment S may be formed by combining information 
from segment S and the previous R - 1 segments, such as 

* * = S 



With values for 



30 



a 



and ^ <s, (//) determined, estimator 70 determines whether there is at least one tone in each of the high and 
low DTMF tone groups which is within ±(1.5% + 2Hz) of the defined DTMF tone frequencies, f t . This is done 
35 by assessing whether the values for 



a 



40 are within the required bounds of f, (such bounds may provide, e.g., that 

ldj^kO.035/,). 

45 The processing described below need be performed only for those tones within the required bounds. If 
there is not at least one tone with 



50 



55 



a<*> 

within the required bounds in both the high and low frequency bands, then processing for the current group 
of segments may cease and the signal slicer 15 may be notified to provide the next segment. 
Once 



is determined, a refined amplitude estimate A {s) (fj) can be determined as follows: 
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AT-1 

L 

k = 0 



Z V k (0)yl*Hfj) + # } *L cM'>(/i) 

Jk = 0 



K-1 

*=0 



V*<0) 



(C.0.2) 



70 



or, more accurately, as follows: 



Z K(0;a^) 
* = 0 1 



(C0.3) 



20 



25 



where 



N- l 

V*(0;a) = X v *(')e 

Often (C.0.2) provides an adequate approximation to A (s) (fj). 
F-Test processor 72 determines F s (fj) as follows: 



i2Ka(r ~/ 0 ) 



30 



(K-l) £ A (/>) 
Fs(fj) - S-(R-l) K.-X 



2 K-\ 

z 

k=Q 



VaCO;*}") 



35 



Z " Z yi'Hfj) ~ A (I \fj)v k {0-a? )\ 

s=S k=0 



(C.0.4) 



This embodiment of the present invention may be augmented as described above with additional 
processing associated with, e.g. a maximum amplitude variation check 52, a maximum phase var.at.on 
check 54, and a twist check 56. 
40 Furthermore, the amplitude and phase variation checks, 52 and 54, may be combined since an estimate 

of the frequency offset in the s' h segment 



a 



45 



is known Consider an illustrative DTMF receiver designed so that R successive segments of N elements 
each have to contain a valid DTMF signal for the receiver to determine that a DTMF signal is present. A 
tone with amplitude B and a frequency near f, can be represented as: 



50 



Be 



.2*/, . t = o,l,2,...N-l 



(C.0.5) 



in the s th segment, where 



55 



fj iS )=fj+CL) 
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is the frequency, and 



a<*> 



10 



15 



20 



25 



is the frequency offset in the s th segment. The initial amplitude of the tone in the s th segment is 
The DTMF receiver produces the amplitude estimates 

which can be phase shifted to yield the quantities £ (s) ,which are estimates of the amplitude B in the s m 
segment: 

B is) = X] s) e- i2nsiN - l)fiU) (C.0.7) 

The true amplitude B can be estimated by an arithmetic average of the B is) : 

B = ± S - { Z X) B {S) = X'-^-^A^e-^- 1 ^ (C.0.8) 
R s=S R s=S 



30 



35 



and the sample variance of the B {s \ 



2 i S-(R-\) 



(C.0.9) 



(^) provides a measure of the spread of the random variables B ls) . The variance ratio 

40 

\B\ 2 

f b = —T (CO. 10) 

45 

is roughly F-distributed with 2 and 2R-2 degrees of freedom. The ratio F B may be compared with a 
threshold, illustratively 2.5. When F B is less than the threshold, the complex amplitude of successive 
segments is determined to vary more than expected for a DTMF tone and the segment may be rejected. 

so D. An Illustrative Embodiment Providing Contemporaneous Fit of Two DTMF Tones 

A further illustrative embodiment of the present invention is presented in Figure 10. This embodiment 
provides for the contemporaneous fit (or estimation) of the, amplitudes of two tones which may form a 
DTMF signal. Unlike the embodiments discussed above, where one or more values for W are chosen to 
55 define limited bands about defined DTMF tone frequencies f h the present embodiment provides a 
frequency band centered about a frequency f c . A value for W is chosen to provide a band about f c which 
includes all eight possible DTMF tone frequencies. In the embodiments previously, described, a model for 
DTMF signal energy (employed in the F-Test) concerned a single frequency component in given frequency 

19 
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bands about frequencies /, In this embodiment, the model for DTMF signal energy includes both frequency 
components representing the low and high tones of a DTMF signal. 

Suppose th'at a frequency band [f„- W.f c * W) is chosen to include all D TMF tone J^ 8 ^ * £ 
exclude call-progress frequencies, f CP (such as those associated w.th d.al tone (350 Hz and 440 Hz)). For 

s example, such a'band may be provided by defining a center frequency of f. 1 500 ^ ^ 

a bandwidth of 2W corresponding to 2000 Hz, so that the band extends from 500 to 2500 Hz 

Ste ^ values for N and W a set of DPSS, , M . K = 0,1,2.K - 1. may be determ.ned to concentrate 
receded signal energy within the band. There are at most \_2NW\oi these sequences (or tapers) which 
esfst teakage. By choosing a wide frequency band, [2NW\ may be a dozen or more. For examp^Mhe 

w duration of signal segment is T = 10 msec, and the sampling penod .s r-1/8000Hz, each segment 
contains N = tU 8 0 samples. If 2W corresponds to 2000 Hz as above, then there are approximate y 
I 2NW\ = K = 20 tapers that have useful leakage properties (i.e., concentration ratios. X,, approximately 
equal to one). Values for the 20 tapers may be determined as described above and stored m memory 26 
alona with the values for tapers as described above for W corresponding 230 Hz. _ 

, 5 L with the embodiments described above, a signa. which may contain a DTMF signal. x(f). ,s provided 
to signal slicer 15, which in turn provides signal segments x«<r). These segments are passed to DFT 
processor 20 which determines eigen-coefficients, 

20 

by performing DFTs as described above using the tapers retrieved from memory 26. 
The eigen-coefficients, 

>4 5) (/A 

are provided to amplitude and frequency estimator 75. In addition, a set of coefficients, 

are determined for later use. Amplitude and frequency estimator 75 determines amplitude estimates, 
as per equation (B.O.2.), and estimates of frequency deviation about f h 



35 



as per equation (C.0.1). With values of 



40 



45 



55 



computed, estimator 75 determines whether there is one and only one tone in each of the low and high 
DTMF tone groups within ±(1 .5% + 2Hz) of the defined DTMF tone frequencies, f„ referred to as 

//and//. 



so res 



pectively. This is done, as in section C. by assessing whether values for 



are within the required bounds of f-„ applying a preliminary F-Test of the form (C.0.4) performing phase ; and 
a^DH ude checks, and a twist check. If other than two single tones (one in each band) are .dent.fied. 
processing of the current segment may cease and the signa, slicer 15 may be notified Otherwise, the 

following processing may be performed. 
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Estimator 75 determines estimates of the frequency of the identified single tones by summing the 
appropriate values of f h i.e., 



fj H and//. 



with associated values of 



w 



&}'\&)'> L and 



15 



fh s) 



= ft + &J' )L 

= f? + *j* )W . 



(£>.0.1) 
(D.0.2) 



20 With frequency estimates determined, estimator 75 determines estimates of amplitude at these frequencies 
in the segment s: 



25 



(D.0.3) 



30 



Ah 



-jr Z >i f) Uc ) 1 V* (fc -fh £) > i V € (fc -ft* ) I 2 - V\ {fc -ft* ) V € (f c -f[ s) ) V\ (/ c -/,$'> ) 



35 



(D.0.4) 



o 



where * represents a complex conjugate and 



40 



K-\ 



V k {f c ~ft s) ) V e (f c -f^) 



45 



- Z Vk(fc-fh s) )K(fc-fl s) )v<(fc-f£ s) )vUfc-f^ oxo.5) 



With values for 

determined, F-Test processor 77 may determine a value for F* s) as follows: 
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yc-2) 
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Z \yi'Hfc)-APv k if e -ftb-AVv k if e -fi'>) 



— . (D O. 6) 



5 



which is F-distributed with 4 and 2K - 4 degrees of freedom, and is a measure of how well the s th segment 
fits the DTMF model. If a segment is associated with a small value of f^, then it probably does not contain 

70 a valid DTMF signal. On the other hand, if the segment is associated with a large value of F< s> , the segment 
would more probably be associated with true DTMF signals. As in sections B and C, the sensitivity of 
(D.0.6) may be enhanced by averaging both numerator and denominator of (D.0.6) over R segments. As 
with the embodiments discussed above, additional tests may be performed to augment this embodiment. 
This matrix-based test, although useful in low signal-to-noise environments, may take more processing 

75 time than other DTMF tests. Therefore, while it may be used generally, it may also be appropriate to limit 
its use to deciding intermediate cases that arise only rarely. The F-Test can be made more discriminating 
by increasing the bandwidth 2 W, as long as this does not cause tones which are a natural part of a 
telephone system, such as call progress tones, to intrude on the frequency band [f c -W t f c + W). If this 
occurs, steps can be taken to make this test useful. For example, the model can be expanded to include 

20 the call progress tones, and the amplitudes of these call progress tones can be determined. However, this 
may increase the size , of the matrix problem that must be solved and lead to increased complexity and 
slower processing. 

E. DTMF Signal Distortion 



DTMF receivers may be confronted with DTMF signals which suffer from distortion, either due to 
problems with a DTMF tone generator, or because of distortion created by a network. For example DTMF 
tones may not have constant frequency, but may suffer from phenomena known as frequency drift or 
chirp. Some telephones may produce DTMF tones which are amplitude modulated, making a single DTMF 
30 signal appear like several repeated DTMF signals. Furthermore, nonlinearities in either a DTMF tone 
generator or a network may result in significant intermodulation components in a DTMF signal. Distortion 
may also take the form of clipping or quantization error. If these sources of distortion are severe, real DTMF 
signals may be rejected by a DTMF receiver at an unacceptable rate. 

In the context of the present invention, these problems may be addressed by relaxing the rejection 
35 thresholds associated with various tests. This, of course, may increase talk-off events. A balance may 
therefore need to be struck in receiver design given the severity of the distortion experienced. 

Alternatively, identified distortion may be introduced into the models of DTMF signals used by the F- 
Test. Tuning the F-Test in this fashion may avoid increasing talk-off events at the expense of complexity. 
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/"Use three fixed dft's (nominal +- 3.5%) per frequency bin */ 
♦include <stdio.h> 
♦include <math.h> 
♦define K 3 
* 5 ♦define PI 3.141592654 
♦define NFREQ 8 
♦define NWIN 3 

♦define USAGE "3dftt3l -nsamples -blinear-3, mu-law sparc*2,mu-law»l / \ 
a"scii~0 -ssJcip ifile dfile ofile" 

♦ define" BUFS2 8192 

♦ define^ FRAT 5 .0 

♦ define' MAXPWR 3.5 
♦define* MINPWR -30.0 
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float f reql (NFREQ } -{ 672 . 0,743. 0,822. 2, 90 8. 1,1166. 7,1289.0,1425. 3, 157 5.8) ; 
float freq [NFREQ] -{ 697 . 0,770. 0,8 52. 0, 94 1.0, 1209. 0,1336.0, 1477. 0, 1633 . 0 } ; 
float f reqh [NFREQ] - { 72 1 . 3 95 , 7 97 . 0 , 8 8 1 . 8 , 97 4 . 0 , 125 1 . 3 , 1 382 . 8 , 152 8 . 7 , 1 6 90 . 2 } ; 
char *tt_table (] [4] - { 

{"1", -2", ** 3" , "A" } , 
{ "4 M , "5", "6", "B" } , 
r7 «,- 8 -, "9 W ,"C~>, 

20 (***-, "0", "I", "D" } 

H ^ * 

main (argc, a rgv) 
int ,argc; 
char " **argv; 

{ ' ? / 

extern int optind; 
extern char *optarg; 
char dtype (8] , ibuf [BUFSZJ ; 
char *ttlogic(); 

double v2, vO [K] , data [ 80 ] , f datar [NFREQ] (K] , f datai (NFREQ] [K] , 

A0r[NFR£Q], AO i [NFREQ] , res [NFREQ] , A^seqr [NFREQ] [NWIN] , 
A_seqi [NFREQ] [NWIN] , res^seq (NFREQ] [NWIN], F (NFREQ] , 
AOlr (NFREQ) , AOli (NFREQ) , AOhr [NFREQ] , AOhi [NFREQ] , 

A_seqlr [NFREQ] [NWIN] , A_seqli [NFREQ] (NWIN) , A_seqhr (NFREQ J [NWIN], 
A_seqhi (NFREQ] (NWIN] , res_seql (NFREQ] (NWIN] , re s_s eqh (NFREQ] [NWIN] , 
35 Fl [NFREQ] , Fh [NFREQ] , real (NFREQ] , reah (NFREQ] , f datalr [NFREQ] [K] , 

f da tali [NFREQ] [K] , f datahr [NFREQ] [K] , fdatahi [NFREQ] ( K] , pow, powdbm, 
nf req [NFREQ] (NWIN] ; 
double taper [80] [K] ; 

int i, n, j, itype, k, rare ad, skip, icnt ; 
FILE *iptr, *optr, *dptr, *fopen () ; 
40 /***************************** set Constants****"************"**** •***/ 
strcpy (dtype, "real m ) ; 
if (argc < 3) 
( 

printf ("%s0, USAGE); 
exit (0) ; 

} 
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while ( ( i-getopt (argc, argv, "n :b: 3 : " ) » EOF) 
switch ( i) 
{ 

case ' n' : 

mread - atoi (optarg) ; 
if (mread > BUFS2) 

1 fprintf (stderr, N Max size - %dO,BUFSZ) 
exit (1) ; 

) 

printf < w %d0, mread); 
break; 
case ' b' : 

itype - atoi (opta rg) ; 
printf r%dO, itype); 
75 break; 

case ' s' : 

skip - atoi (optarg) ; 
printf ("%d0, skip) ; 
break; 

case * ~* ' • 

fprintf (stderr, "ERROR: Illegal optionO) ; 
exit (1) ; 

/*******i***********open files calculate statistics 
iptr - fopen (argv [opt.ind++] , "r") ; 
dptr - fopen (argv [optind++3 # "r") ; 
optr - fopen (argv (opt ind++] , "a") 
if (itype— 2 II itype— 1) 
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if (f seek (iptr , skip, 0) ) 

( 

perror ( w f seek") ; 
exit ID-; 

} 

) 

else if (itype— 3) 

if (f seeMiptr, skip*sizeof (short) , 0) ) 
{ 

perror (*f seek") ; 
exit(l); 

} 



else if(skip!-0) 

{ printf ("Error : Ascii Mode with non-zero skipO) ; 
45 exit (1) ; 

} 

/-Collect DPSSs */ 
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for (i-0; fscanf (dptr, #, %lf %lf %lf %* 1 f ,& tape r t i 1 C 0 \, & tape r [ i ][ 1 ] , 

&taper [i] (2) ) !« ECF;i++); 
/* Compute taper constants "/ 
v2-0.0; 

for <k-0;k<K;k++) 

{ 

vO[k]-0.0; 

' f or ( i a 0 ; i<mread; i + + ) 
( 

vO [k] +-taper [ij (k] ; 



} 

v2+-v0 [kl *v0 {k] 



} 



icnt-0; 

for (i-0; KKFREQ; i + + ) 

t 

for ( j-0; j<NWIN; j++> 
{ 

A_seqlr f i] t j] -0 .0; 
A_seqhr[i] (j')-0.0; 
1 A_seqli[i] [j]-0.0; 

A_seqhi[i] [ j]-0 .0; 
v res_seql { i] ( j ] -0 . 0 ; 

res_seqh[i] [jJ-0.0; 
\ A_seqi[i] ( j]-0 .0; 
A_seqr (i] [ j]-0 .0; 
res_seq ( i] { j ] -0 . 0 ; 

} 

} 

30 while (read__blk (mread, itype, data/ iptr) ) 

{ 

icnt++; 

if (itype- -0) 

{ 

for ( i-0 ; i<mread; i++) 
( 

data (i]-4004. 19*data [i] ; /*Convert to standard dbm 

} 

} 

/♦Energy Detector*/ 
40 pow-0 . 0; 

for (i-0; Kmread; i++) 
{ 

pow+-data [i] *data ( i] ; 

) 

pow/- (double) (mread) ; 

powdbm-10 . 0*logl0 (pow) -20 . 0*logl0 (4004 . 19) ; 
if (powdbm < MINPWR) - 
( 
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fprintf (opt r, "%s %ci0 , "Non-ene rgy in this window" , icr.t ) ; 

} 

else 
{ 

if (powdbm > KAXPWR) 
{ 

fprintf (optr, tt %s %d0 , "Possible Clipping" , lent ) ; 

) 

/* Compute Eigencoef f icients */ 

for (n-0 ; n<NFR£Q; n + +) 

( 

for (k-0;k<K;k++) 
( 

fdatalr [n] [k]-0 .0; 
fdataliln] [k]-0.0; 
fdatar [n] [k]«0.0; 
fdatai(n) (k]-0.0; 
fdatahr (n] [k]-0.0; 
fdatahifn] [k]-0.0; 
for ( i«0 ; Kmread; i++) 
( 

fdatalr [n] (k] +- cos<2,00 *PI *i* f reql [n] / 8000 . 0 ) 

*data[i] *taper(i] [k]; 
fdataliln] Ik]+- -sin<2.00 *PI *i* f reql [n] / 8000 . 0 ) 

*data[i] *taper{i] [k] ; 
fdatar [n] [k]+- cos<2.00 *PI *i* f req [n] / 8000 . 0 ) . 

*datafi] *taper [i] [k] ; 
fdatai[n] [k]+- -sin<2.00 *PI *i* f req(n]7 8000 . 0 ) 

*data[i] *taper ( i] (k] ; 
fdatahr {n] [k]+- cos(2.00 *PI *i* f reqh [n ] / 8000 . 0 ) 

*data[i] *taper(i] [k] ; 
fdatahifn] [k}+- -sin(2.00 *PI *i* f reqh [n] / 8000 . 0 ) 

*data[i]. *taper[i] [k] ; 

} 

} 

) 

/* Compute first guess of amp AO (f] and residuals */ 
initguess (fdatar, f datai, vO , v2 , AOr, AOi, res); 
initguess (fdatalr, fdatali, vO, v2, AOlr, AOli, resl) ; 
initguess (fdatahr, f datahi, v0 , v2 , AOhr, AOhi, resh) ; 
/ *Comput« F-values*/ 

f atest <v2, A — seqr, A_ seqi, res_seq, F) ; 
fatest (v2,A_ seqhr, A_seqhi, res_seqh,Fh) ; 
f atest (v2,A_ seqlr, A_seqli, res_seql,Fl) ; 

/♦Decide if DTMF or not based on A_seq, res_seq, nfreq,FA*/ 
fprintf (optr,"%3 %d0, tt logic ( A_seqhr, A_seqhi, res_seqh, A_secl 
,A_seqli, res^seql, A_seqr, A_seqi, res_seq, Fh, Fl , F, v2 ) , ic 
f flush (optr) ; 
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) 

read_blk (size, mode / bu*f , f p) 
mt aire, mode ; 
FILE *fp; 
double bu f [ ) ; 

{ 

mt i; 
float x; 

char ibuffBUFSZ]; 
short sbuf ( BUFSZ ] ; 

switch (mode) 
{ 

case C : 

for ( i-0 ; i<size ; i++) 
{ 

if (facanf (fp, "Sslf &buf ( i] ) — EOF) 
( 

return (0) ; 

} 

} 

return (size) ; 

case 1: 

if <f read (ibuf , sizeof (char) /3ize,fp) r-size) 
{ 

return (0) ; 

} 

else 
{ 

for (i-0; i<aize;i++) 
{ 

buf [i] - (double) mutol (ibuf [i] ) ; 

} 

return (size) ; 

} 

case 2: 

if <fread(ibuf ,sizeof (char) ,size, fp) !«size) 
return (0) ; 

> 

else 

< 

for (i-0; i<size; i++) 
< 

buf [i]- (double) mutol (-ibuf [i] ) ; 

) 

return (size) ; 

) 
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case 3: 

if <f readObuf , sizeof < abort: ) , s i ze, i p) !-aize) 

. ( 

return (0) ; 

} 

else 
{ 

f or ( i-0 ; Ksize ; i + + ) 
( 

buf (i3- (double) (sbuf [x] ) ; 

} 

return (size) ; 

} 



) 



initguess (ar f ai, v, v2,xr,xi, res) 

double art][K) # ai [ ] [K] , v[ ] , v2,xr [] , xi [] , res El ; 

( 

double tmpr,tmpi; 
int i,j; 



for (i-0; i<NFR£Q; i + + > 
i 

tmpr-0 . 0 ; 
tmpi-0 . 0 ; 

25 for ( j-0; j<K; j++) 

( 

tmpr+-ar[i] [ j J *v( j] ; 
tmpi+-ai{i] [j]*v[j]; 



xr [i] «tmpr/v2; 
xi [i] -tmpi/v2; 
tmpr-0.0; 
tmpi-0. 0; 
res [i] -0 .0; 
for ( j-0; j<K; j++) 
{ 

tmpr-ar [i] [ j] -xr ( i] *v t j ] 
tmpi«ai[i] [ j 3 -xi ( i] *v [ j ] 
tmpr *-tmpr ; 
tmpi*-tmpi; 
res [ i ] +-tmpr +tmpi ; 

40 ) 
) 

return; 



fate3t<v2,A_seqr,A_3eqi,re3_seq,F) r , P »« JTM , PM 

double v2,A_seqr(] (NWIN) , A_seqi [ 3 [NWIN] , res_seq [ 3 [NWIN3,F(3 
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double tmpn,tmpd; 
int i, j ; 

for (i-0; i<NFR£Q; i++) 

{ 

cmpn-tmpd-0 . 0 ; 
f or ( j-0 ; j<KWIN; 
{ 

tmpd+- re s_aeq ( i ] ( j ] ; 

tmpn + -A_seqr ( i ] [ j ] *A_seqr ( i ] [ j ] +A_aeqi ( i ] [ j ] *A_aeqi ( i ] [ j J ; 

} 

F[i]- (double) (-1+K) *tmpn *v2 /tmpd; 

} 

return; 



15 char *tt logic ( A_ aeqhr, A^aeqhi, rea — seqh, A_seqlr , A_ seqli, rea_seql, A — seqr , A_seqi, 

res_seq, Fh,Fl, F, v2) 
double A_seqhr( 3 (NWIN] , A_aeqhi [ ] (NWIN] , rea^seqh ( ] (NWIN] , 

A_aeqlr [ ] [NWIN] , A_seqli ( ) (NWIN] , res_seql ( ] (NWIN) , Fh ( ] , Fl ( J , v2 , 
A_3eqr[] (NWIN] , A_seqi ( ] (NWIN] , res_aeq( ] [NWXN],F(]; 
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int i, cntl, cnth, iL, iH, cnt, cntf 1, cntf h; 

double m3lr,mali, aal2, a3h2,m3hr,m3hi,3lr (NWIN] , sli (NWIN] , 

ahr (NWIN] , shi (NWIN] , na, vl, vh, val r vah, FAL, FAH, mreal, mreah,. 
FRL, FRH, pwl (NWIN] , pwh (NWIN] , tw, mash, masl, Fconat , ffvl, fvh, 
FPL, FPH, £1/ £h/ Air (NWIN] , All (NWIN] , Ahr (NWIN] / Ahi (NWIN] ,mAl,mAh, 
m3qAl,m3qAh,absAl (NWIN) , aba Ah (NWIN] , Nr# Ni, Dr, Di, prat 1/ prath, 
resl (NWIN) , resh (NWIN] , vrl, vrh, amph, ampl, amp; 



/*F-teat*/ 
cnt 1-0; 

for <i-0;i<NFREQ/2;i++> 
{ 

if ( F(i) > FRAT) 
( 

cntl++; t 
iL-i; 

> 

35 } 

cnth-0; 

f or < i-NFREQ/2 ; i<NFR£Q; i++ ) 

( 

if (F(i) > FRAT) 
{ 

4 <> cnth++; 

iH-i; 

.) 

) 

/*Only one low 4 high frequency ?*/ 
if <cntl!-l I I cnth!-l) 
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return ("Failed F-test") ; 

el3e /* Passed F-test •/ 

( 

/*Max F value in the midband*/ 
for (i-0; KNFREQ/2; i++) 
{ 

if(il-iL) 
{ 

if((Fh(i) > FRAT) ||<Fl[i] >FRAT) ) 

return ("Failed Sideband F-Teat") ; 

) 

) 

15 ) 

for (i-NFREQ/2; KNFREO; i + + > 
{ 

if (i!-iH> 

( if((Fh[i] > FRAT) M(Fl(i] >FRAT) ) 

20 , 

returnCFailed Sideband F-Test") ; 

) 

) 

25 /*Midband Test for amplitudes*/ 

amp-amp l-amph-0 . 0 ; 
cntf l-cntfh-0; 
f or (i-0 ; KNWIN; 

( ampl-A seqlr[iL] [ i] *A_seqlr [ iL] [i] +A_3eqli[iL] [i] *A_seqli [iL] [ 
30 amph-A~seqhr[iL] [i] *A_seqhr [ iL] [i] +A_seqhi[xL] [i] •A_seqhi UL] i 

amp-A_7eqr[iL] [i] *A_seqr { iL] [i] +A_seqi[iL] [ i] * A_seqx C iL] [i] ; 
if ( ( (amp/ampl) <1 . 00) II ( (amp/amph) <1.00)) 
{ 

cntfl++; 

35 > 
} 

for ( i-0 ; KNWIN; i++) 

1 ampl-A_seqlr[iH] [ i] *A_seqlr [ iH] [i] 

+A_seqli[iH] [i] *A_seqli[iH] [i] ; 
amph-A__seqhr (iH] [ i] *A_seqhr [ iH] {i] 

+A_seqhi(iH] ( i] *A_seqhi C iH] [i]; 
amp-A_seqr [iH] { i ] *A_seqr [ iH ] (i] 

+A__seqi(iH] [i] *A_aeqi(iH] [i] ; 
if (((amp/ampl) <1. 00) II ((amp/amph) <1.00)) 
{ 

cntf h++; 

} 
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) 

if<cntfl >2 II cr.tfh > 2) /"At least one segment has hi n : ^- 

( 

return ( Tailed Mid-Band Amplitude Test"); 

I 

lor ( i-0; i<NWIN; i + +) 
( 

Air [i]-A_3eqr[iL] [i] ; 
Ali[i)«A_seqi[iL] [i] ; 
fl-freq(iL); 
res 1 [ i} «res_seq [ iL] (i]; 
Ahr (i)«A_seqr (iH3 [i] ; 
Ahi [i)«A_seqi (iH) [i] ; 
fh-f reqli] ; 

resh ( i) -res_seqh [ iH] (i]; 

is ) 

mAl-mAh-msqAl~msqAh-0 . 0 ; 

for (i-0; KNWIN; i + + ) 
( 

absAl[i)-sqrt (Air (i) *Alr (ij +Ali(i] *Ali(i)) ; 
abs Ah ( i ] -sqrt < Ahr ( i ] * Ah r ( i J +Ahi [ i ) *Ahi ( i] ) ; 
20 mAl+-absAl(i) ; 

mAh+-absAh ( i] ; 
msqAl+-absAl {i] *absAl f i] ; 
msqAh+-abs Ah { i ] * abs Ah [i] ; 

} 

na- (double) (NWIN) ; 
mAl/-na; 
mAh/«na; 

val-(msqAl - na*mAl*mAl) / <na-l . 00) ; 
van- (msqAh - na*mAh*mAh) / (na-1 . 00) ; 
/*F-test for amplitude variance*/ 
FAL-2 .00*mAl*mAl /val; 
FAH-2 .00*mAh*mAh /vah; 
if (FAL <FRAT I I FAH <FRAT> 
( 

return ( "Excessive amplitude variation"); 

>* 

35 else 
{ 

/♦Frequency and phase jump test*/ 

Nr-Alr [0 ] * Alr.[ 2 ] -Ali [ 0 ) *AXi C 2 ] -Air ( 1 ] *Alr [ 1 ] +Ali [ 1 ] * Ali [ 1 J ; 
Ni-Ali ( 0 3 * Air [ 2 ) +Alr { 0 ] *Ali ( 2 ] -2 *Alr ( 1 ) * Ali ( 1 ] ; 
Or-Alr(l] *AlrC2]-Ali(l] * Ali 12] ; 
40 Di-Alr(l]*AliC2]+Ali(l)*Alrt2] ; 

pratl-(Nr*Nr+Ni*Ni) / (Dr *Dr+Di*Di) ; 

Nr-Ahr[0]*Ahr{2]-Ahi[0] *Ahi [2 ] -Ahr [ 1 3 *Ahr [ 1 3 +Ahi [1 ] *Ahi [ 1 ] ; 
Ni-Ahi(0 3 *Ahr (2]+Ahr [03 *Ahi (2 ] -2* Ahr [ 1 3 *Ahi [1] ; 
45 Dr-Ahr(l] *Ahr (2 ) -Ahi [ 1 3 *Ahi(2] ; 

Di-Ahr(l] *Ahi[2)+Ahi[l] *Ahr [2] ; 
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prath-<Nr«Nr+Ni«Ni)/<Dr-Dr*Di-Di> ; 

/•Phase jump test'/ 

if <<pratl > 0.4) l I (prath >0.4>> 

1 return ("Excessive phase jump"); 

5 

} 

else 
{ 

/*twisr test*/ 
for(i-0;i<NWIN;i++) 

owl [ i] -absAl ( i] *absAi [ i] ; 
pwh t i ] *abs Ah ( i ] * abs Ah U 1 ; 
tw-10.0*logl0 (pwllil /pwhUl ) - 
if ((tw<-6.0) I I <tw>8.0>) 

^ ( re turn ("Excessive Twist") 

) 

) 



/* Valid DTMF Signal */ 
return ( tt_t able [iL] (iH-4]> ; 
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Claims 
1. 



system channel, the method comprising the steps of: 

applying a plurality of tapers to a segment of a communicated signal; 

communication system element. 

The method of claim 1 further comprising the step of receiving a communicated signal and producing 

one or more segments thereof. 

The method of claim 2 wherein the step of producing one or more segments comprises the step of 
producing overlapping segments. 

The method of claim 1 wherein the control signal information comprises DTMF signal information. 
The method of claim 1 wherein the contro! signal information comprises call progress signal informa- 



40 2. 



4. 

5. 



50 



tion. 

6. 



55 



The method of claim 1 wherein the control signal information comprises supervisory audio tone signal 

information. 

7 The method of claim 1 wherein the step of applying a plurality of tapers comprises the step of applying 
a taper which comprises a discrete prolate spheroidal sequence. 

8 The method of claim 1 wherein the step of applying a plurality of tapers comprises the step of forming 
a product of a taper value and an input signal segment value. 
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9. The method of claim 1 wherein the step of applying a plurality of tapers comprises the step of applying 
a taper value to a possible communicated signal value. 

10. The method of claim 9 further comprising the step of storing in memory information which comprises a 
5 taper value applied to a possible communicated signal value. 

11. The method of claim 1 wherein the step of performing a transform comprises the step of performing a 
Fourier transform. 

.10 12. The method of claim 1 wherein the model of a control signal comprises a sinusoidal model. 

13. The method of claim 12 wherein the sinusoidal model comprises an amplitude modulated sinusoidal 
model. 

75 14. The method of claim 12 wherein the sinusoidal model comprises a frequency modulated sinusoidal 
model. 

15. The method of claim 12 wherein the sinusoidal model comprises a first sinusoidal model and a second 
^ i sinusoidal model. 

~" 20 

16. The method of claim 15 wherein the first sinusoidal model comprises a first control signal information 
frequency and the second, sinusoidal model comprises a second control signal information frequency. 

17. The method of claim 16 wherein the step of determining a similarity score comprises the step of 
25 estimating amplitude information based on a transformed tapered segment and wherein each of the first 

and second sinusoidal models comprises estimated amplitude information. 

18. The method of claim 17 wherein the step of determining a similarity score comprises the step of 
forming a ratio, the numerator of the ratio comprising an energy distribution based on the first and 

30 second sinusoidal models, and the denominator of the ratio comprising an energy distribution based on 

a difference between a transformed tapered segment and a sum of the first and second sinusoidal 
models. 

19. The method of claim 1 wherein the model of a control signal comprises expected control signal 
35 information distortion. 

20. The method of claim 1 wherein the step of determining a similarity score comprises the step of 
estimating amplitude information based on a transformed tapered segment and wherein a model of a 
control signal comprises estimated amplitude information. 

40 

21. The method of claim 1 wherein the step of determining a similarity score comprises the step of 
determining a measure of similarity between energy distributions. 

22. The method of claim 21 wherein the step of determining a measure of similarity between energy 
45 distributions comprises the step of forming a ratio, the numerator of the ratio comprising an energy 

distribution based on a control signal model, and the denominator of the ratio comprising, an energy 
distribution based on a difference between a transformed tapered segment and the control signal 
model. 

50 23. The method of claim 1 wherein the step of determining a similarity score comprises the step of: 

determining an estimate of frequency deviation between an expected control signal information 
frequency and a frequency component of a transformed tapered segment; and 
determining whether the estimate of frequency deviation exceeds a threshold. 

55 24. The method of claim 1 wherein the step of identifying control signal information comprises the step of 
determining whether a similarity score exceeds a threshold. 
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25. The method of claim 1 wherein the step of identifying control signal information comprises the step of 
determining whether a frequency component of a transformed tapered segment deviates from an 
expected control signal information frequency by an amount which exceeds a threshold. 

» 26. The method of claim 25 wherein the step of determining whether a frequency component of a 
transformed tapered segment deviates from an expected control signal information frequency com- 
prises the steps of: 

determining amplitude estimates of the transformed tapered segment at a plurality of frequencies 
including the expected control signal information frequency; and 
o determining whether the amplitude at the expected control signal information frequency exceeds 

the other determined amplitude estimates. 

27. The method of claim 1 wherein the step of identifying control signal information comprises the steps of: 

determining an amplitude variation parameter based on estimates of the amplitude of a frequency 
5 component present in each of a plurality of segments; and 

determining when the amplitude variation parameter exceeds a threshold. 

28. The method of claim 1 wherein the step of identifying control signal information comprises the steps of: 

determining a phase variation parameter based on estimates of the amplitude of a frequency 
»o component present in each of a plurality of segments; and 

determining when the phase variation parameter, exceeds a threshold. 

29. The method of claim 1 wherein the step of identifying control signal information comprises the steps of: 

determining an amplitude and phase variation parameter based on estimates of the amplitude and 
?5 phase of a frequency component present in each of a plurality of segments; and 

determining when the amplitude and phase variation parameter exceeds a threshold. 

30. The method of claim 1 wherein the step of identifying control signal information comprises the steps of: 

determining a ratio of amplitudes of a transformed tapered segment; the numerator of the ratio 
30 comprising an amplitude evaluated at one frequency of expected control signal information and the 

denominator of the ratio comprising an amplitude evaluated at another frequency of the expected 
control signal information; and 

determining whether the ratio is within a predetermined range of values. 

35 31. An apparatus for detecting control signal information for a communication system element, the control 
signal information included in one or more segments of a signal communicated over a communication 
system channel, the apparatus comprising: 

means for applying a plurality of tapers to a segment of a communicated signal; 
means for performing a transform of a tapered segment; 
40 means for determining a similarity score based on a transform of a tapered segment and a model 

of a control signal; 

means for identifying control signal information based on a similarity score; and 
means for providing an indicator signal representative of the identified control signal information to 
the communication system element. 

45 

32. The apparatus of claim 31 further comprising a means for receiving a communicated signal and 
producing one or more segments thereof. 

33. The apparatus of claim 32 wherein the means for producing one or more segments comprises means 
so for producing overlapping segments. 

' 34. The apparatus of claim 31 wherein the control signal information comprises DTMF signal information. 

35. The apparatus of claim 31 wherein the control signal information comprises call progress signal 
55 information. 

36. The apparatus of claim 31 wherein the control signal information comprises supervisory audio tone 
signal information. 
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37. The apparatus of claim 31 wherein the means for applying a plurality of tapers comprises a means for 
applying a taper which comprises a discrete prolate spheroidal sequence. 

38. The apparatus of claim 31 wherein the means for applying a plurality of tapers comprises a means for 
5 forming a product of a taper value and an input signal segment value. 

39. The apparatus of claim 31 wherein the means for applying a plurality of tapers comprises a means for 
applying a taper value to a possible communicated signal value. 

io 40. The-apparatus of claim 39 further comprising a memory for storing information which comprises a taper 
value applied to a possible communicated signal value. 

41. The apparatus of claim 31 wherein the means for performing a transform comprises a means for 
performing a Fourier transform. 



75 



42. The apparatus of claim 31 wherein the control signal model comprises a sinusoidal model. 



43. The apparatus of claim 42 wherein the sinusoidal model comprises an amplitude modulated sinusoidal 
^ ) model. 
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44. The^apparatus of claim 42 wherein the sinusoidal model comprises a frequency modulated sinusoidal 
model. 

45. The ^apparatus of claim 42 wherein the sinusoidal model comprises a first sinusoidal model and a 
25 second sinusoidal model. 

46. The apparatus of claim 45 wherein the first sinusoidal model comprises a first control signal information 
frequency and the second sinusoidal model comprises a second control signal information frequency. 

30 47. The apparatus of claim 46 wherein the means for determining a similarity score comprises means for 
estimating amplitude information based on a transformed tapered segment and wherein each of the first 
and second sinusoidal models comprises estimated amplitude information. 

48. The apparatus of claim 47 wherein the means for determining a similarity score comprises means for 
35 forming a ratio, the numerator of the ratio comprising an energy distribution based on the first and 

second sinusoidal models, and the denominator of the ratio comprising an energy distribution based on 
^ a difference between a transformed tapered segment and a sum of the first and second sinusoidal 

r^) models. 
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49. The apparatus of claim 31 wherein the control signal model comprises expected control signal 
information distortion. 



50. The apparatus" of claim 31 wherein the means for determining a similarity score comprises a means for 
estimating amplitude information based on a transformed tapered segment and wherein a control signal 

45 model comprises estimated amplitude information. 

51. The apparatus of claim 31 wherein the means for determining a similarity score comprises means for 
determining a measure of similarity between energy distributions. 

50 52. The apparatus of claim 51 wherein the means for determining a measure of similarity between energy 
distributions comprises means for forming a ratio, the numerator of the ratio comprising an energy 
distribution based on a control signal model, and the denominator of the ratio comprising an energy 
distribution based on a difference between a transformed tapered segment and the control signal 
model. 
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53. The apparatus of claim 31 wherein the means for determining a similarity score comprises: 

means for determining an estimate of frequency deviation between an expected control signal 
information frequency and a frequency component of a transformed tapered segment; and 
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means for determining whether the estimate of frequency deviation exceeds a threshold. 

54. The apparatus of claim 31 wherein the means for identifying control signal information comprises 
means for determining whether a similarity score exceeds a threshold. 

55. The apparatus of claim 31 wherein the means for identifying control signal information comprises 
means for determining whether a frequency component of a transformed tapered segment deviates 
from an expected control signal information frequency by an amount which exceeds a threshold. 

w 56. The apparatus of claim 55 wherein the means for determining .whether a frequency component of a 
transformed tapered segment deviates from an expected control signal information frequency com- 
prises: 

means for determining amplitude estimates of the transformed tapered segment at a plurality of 
frequencies including the expected control signal information frequency; and 
75 means for determining whether the amplitude at the expected control signal information frequency 

exceeds the other determined amplitude estimates. 

57. The apparatus of claim 31 wherein the means for identifying control signal information comprises: 

means for determining an amplitude variation parameter based on estimates of the amplitude of a 
20 frequency component present in each of a plurality of segments; and 

means for determining when the amplitude variation parameter exceeds a threshold. 

58. The apparatus of claim 31 wherein the means for identifying control signal information comprises: 

means for determining a phase variation parameter based on estimates of the amplitude of a 
25 frequency component present in each of a plurality of segments; and 

means for determining when the phase variation parameter exceeds a threshold. 

59. The apparatus of claim 31 wherein the means for identifying control signal information comprises: 

. means for determining an amplitude and phase variation parameter based on estimates of the 
30 amplitude and phase of a frequency component present in each of a plurality of segments; and 

means for determining when the amplitude and phase variation parameter exceeds a threshold. 

60. The apparatus of claim 31 wherein the means for identifying control signal information comprises: 

means for determining a ratio of amplitudes of a transformed tapered segment, the numerator of 
35 the ratio comprising an amplitude evaluated at one frequency of expected control signal information 

and the denominator of the ratio comprising an amplitude evaluated at another frequency of the 
expected control signal information; and 

means for determining whether the ratio is within a predetermined range of values. 

40 
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